Large-g asymptotics of the random bond Potts model 
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ABSTRACT 

We numerically examine the large-q asymptotics of the g-state random bond 
Potts model. Special attention is paid to the parametrisation of the critical line, 
which is determined by combining the loop representation of the transfer ma- 
trix with Zamolodchikov's c-theorem. Asymptotically the central charge seems 

o\ ; 

Q\ \ to behave like c{q) = 2log2('z) + 0{1). Very accurate values of the bulk mag- 

C^ ' netic exponent xi are then extracted by performing Monte Carlo simulations di- 

rectly at the critical point. As g ^ oo, these seem to tend to a non-trivial limit, 

Ch ; xi ^ 0.192 ± 0.002. 

O 

O 
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1 Introduction 

Recently the two-dimensional g-state random bond Potts model with q > 4 has attracted 
considerable interest, because it serves as a paradigm for examining the effect of quenched 
randomness []I| on a first-order phase transition 0. Since in this case the randomness 
couples to the local energy density, a theorem by Aizenman and Wehr 0, along with 
related analytical work 0, ||, suggests that the transition should become continuous, 
as has indeed been verified by subsequent numerical studies [^, |^, ||, ^ [1^, |1T], |T2|, |13 



Unfortunately, analytical results have been scarce, except in the limit g ^ oo where 
properties of a particular tricritical point were related to those of the zero-temperature 
fixed point of the random field Ising model in d = 2 + e dimensions p[. From the con- 
jectured phase diagram [B| it is however known that this fixed point is not the analytical 



continuation of the line of random fixed points found for finite q > 2 |TJ, |T5| . Namely 
the latter (henceforth referred to as the q ^ oo limit of the model) is rather believed to 
be associated with a subtle percolation-like limit P], the exact properties of which have 
not yet been fully elucidated. 

In the present publication we seek to gain further knowledge of this q ^ oo limit 
by producing numerical results along the afore-mentioned line of critical fixed points for 
very large values of q. Since cross-over effects to the pure and percolative limits of the 
model have been shown to be important p, |I^, special attention must be paid to the 



parametrisation of the critical line. Generalising a recently developed transfer matrix 



technique |jT6[, in which the Potts model is treated through its loop representation [|T7 
we were able to explicitly trace out this line, and as a by-product obtain very precise 
values of the central charge. Based on our numerical results for the q = 8^ state model 
with k = 1,2, ... ,6 we find compelling evidence that 

c{q) = hog,{q) + 0{l). (1.1) 

Although this behaviour of the central charge is reminiscent of the Ising-like features 
of the tricritical fixed point discussed above, we shall soon see that from the point of 
view of the magnetic exponent the q ^ oo limit is most definitely not in the Ising 
universality class. Note also that our precision allows us, for the first time, to convincingly 
distinguish the numerically computed central charge from its analytically known value 

1 



in the percolation limit 0. 

With the numerically obtained parametrisation of the critical disorder strength at 
hand we then proceed to measure the corresponding magnetic bulk scaling dimension xi 
as a function of q. The most suitable technique is here that of conventional Monte Carlo 
simulations. Our results lend credibility to the belief [O that Xi(g) saturates as g -^ oo. 



Based on results for the q = 8'' state model with k = 1, 2, 3 we propose the limiting value 

xi(g) ^0.192 ±0.002 for q ^ oo, (1.2) 

in agreement with the one reported in Ref . [|l^ . The fact that Eq. (|1.2| ) does not coincide 
with any known scaling dimension of standard percolation is remarkable, and calls for 
further analytical investigations of the q ^ oo limit. 

After explaining the loop model transfer matrices in Section ^ we state our results 
for the critical line and the central charge in Section ^. The Monte Carlo method and 
the resulting values of the magnetic scaling dimension are presented in Section ^ and we 
conclude with a discussion. 

2 Loop model transfer matrices 

The partition function of the random bond Potts model can be written as 

Z = Y,Yle^"^"-"': (2.1) 

W} {ij) 

where the summation is over the q discrete values of each spin and the product runs 
over all nearest-neighbour bonds on the square lattice. The Kij are the reduced coupling 
constants, which for the moment may be drawn from an arbitrary distribution. By the 
standard Kasteleyn-Fortuin transformation |18|, Eq. ( p.lD can be recast as a random 
cluster model 

Z = T.f'^'^ n (e^'^-l), (2.2) 

{G} {ij)(^g 

where ^ is a bond percolation graph with C{Q) independent clusters. Note that q now 
enters only as a (continuous) parameter, and since the non-locality of the clusters does 



not obstruct the construction of a transfer matrix |jl9[ the interesting regime of g ^ 4 



becomes readily accessible, provided that one can take into account the randomness in 
the couplings 0. 

In an analogous fashion we can adapt the even more efficient loop model representation 
||T6| to the random case. Indeed, trading the clusters for their surrounding loops on the 
medial lattice [|r^, Eq. ( [^.2[ ) is turned into 

^ = ,A^/2^^L(g)/2 J^ f^^), (2.3) 

{g} (ij)eg V V^ / 

where N is the total number of spins, and configuration Q encompasses L{Q) loops. The 
strip width L is measured in terms of the number of 'dangling' loop segments, and must 



be even by definition of the medial lattice |16 



A pleasant feature of the random bond Potts model is that the critical temperature 
is known exactly by self-duality pOf. Employing for simplicity the bimodal distribution 



P{Ki^) = i [6{K,, - K,) + 6{Ki, - K^)] , (2.4) 

and choosing the parametrisation Sij = (e^'^ — I)/ ^/q, the self-duality criterion takes the 
simple form 

siS2 = 1. (2.5) 

To fully identify the critical point the only free parameter is then the strength of the 
disorder, which can be measured in terms oi R = K1/K2 > 1 or s = si > 1. 

3 Central charge 



In Ref. |T^ we showed that Zamolodchikov's c-theorem pl| is a powerful tool for nu- 
merically identifying the fixed points of a pure system. The idea is simple: From the 
leading eigenvalue of the transfer matrix, specific free energies fo{L) can be computed 
as a function of the strip width L. Effective central charges c{L) are then obtained by 
fitting data for two consecutive strip widths according to p2| 



/o(^) = /o(oo)-^ + ---. (3.1) 

By tuning the free parameter s of the system, local extrema c{L, s*(L)) are sought for, 
and finally the fixed point is identified by extrapolation: s* = s*(L — > 00). 



In principle this strategy can also be employed for a disordered system, provided that 
error bars are carefully kept under control. Now fo{L) is related to the largest Lyapunov 



exponent of a product of Af ^ oo random transfer matrices p3|, §, and its statistical 
error vanishes as M~^/^ by the central limit theorem. Thus, for large enough M any 
desired precision on fo{L) can be achieved. 

An important observation is that for larger and larger L, the c{L) found from Eq. (|3.1| ) 
become increasingly sensitive to errors in fo{L). Therefore M must be chosen in accor- 
dance with the largest strip width Lmax used in the simulations. For the system at hand 
we found that 4 significant digits in c{L) were needed for a reasonable precise identifica- 
tion of s^:{L), and with L^ax = 12 this in turn implies that the fo{L) must be determined 
with 6 significant digits. We were thus led to choose M = 10^ for q = 8, and M = 10^ 
for larger values of g.0 

Data collection was done by dividing the strip into M/l patches of length / = 10^ lat- 
tice spacings, and for each patch the couplings were randomly generated from a canonical 
ensemble, i.e., the distribution (|2.4|) was restricted to produce an equal number of strong 
and weak bonds. 

In the right part of Table ^ we show the resulting two-point fits ( |3.1| ) in the q = 8 
state model, as a function of s. The left part of the table provides analogous three- 
point fits, obtained by including a non-universal 1/L^ correction in Eq. ( |3.1| ). In all 
cases the error bars are believed to affect only the last digit reported. The two-point 
fits give clear evidence of a maximum in the central charge, and we estimate its location 
as s* = 6.5 ± 2.0. The corresponding central charge is estimated from the three-point 
fits, as these are known to converge faster in the L —>■ oo limit [Q, and we arrive at 
c = 1.530 ± 0.001. To appreciate the precision of this result, we mention that the 
numerical values of c(g = 8) first reported were 1.50 ± 0.05 |0 and 1.517 ± 0.025 0. 

Table ^ summarises our results for other values of q. Two remarkable features are 
apparent. First, s* oc q"^ is well fitted by a power law with w = 0.31 ± 0.02. This gives 
valuable information on how the g ^ oo limit of the model is approached, and implies 
"^Incidentally, improving our results to Lmax = 14 would require augmenting M by at least a factor 
of 100 (apart from the increased size of the transfer matrices), and since several months of computations 
were spent on the present project this hardly seems possible in a foreseeable future. 



s 


c(4,8) 


c(6,10) 


c(8, 12) 


c(4,6) 


c(6,8) 


c(8,10) 


c(10,12) 


3 


1.495 


1.500 


1.500 


1.4101 


1.4544 


1.4731 


1.4821 


4 


1.512 


1.517 


1.516 


1.4157 


1.4657 


1.4868 


1.4967 


5 


1.519 


1.525 


1.523 


1.4152 


1.4690 


1.4918 


1.5025 


6 


1.521 


1.528 


1.527 


1.4116 


1.4683 


1.4927 


1.5044 


7 


1.520 


1.529 


1.529 


1.4067 


1.4656 


1.4915 


1.5041 


8 


1.518 


1.528 


1.529 


1.4013 


1.4619 


1.4890 


1.5026 


9 


1.509 


1.530 


1.528 


1.3972 


1.4552 


1.4860 


1.5004 


10 


1.511 


1.525 


1.527 


1.3908 


1.4534 


1.4826 


1.4977 


11 


1.501 


1.526 


1.526 


1.3873 


1.4465 


1.4791 


1.4949 


12 


1.504 


1.519 


1.524 


1.3816 


1.4451 


1.4756 


1.4919 



Table 1: Effective central charge of the q = 8 state model, as a function of disorder 
strength s. Two- and three-point fits to Eq. ( p.l| ) are labelled as c{L, L+2) and c{L, L+4) 
respectively. 

that the ratio of the coupling constants R = K1/K2 = log(l + s^/q) / \og{l + ^/q/s) is a 
non-monotonic function of q that tends to the finite limiting value jz^ = 4.3 ± 0.6 as 
g ^ 00. We shall discuss this finding further in Section ^ Second, the central charge 
seems to fulfil the relation (|TTT|) as stated in the Introduction. 

4 Magnetic scaling dimension 



In this section we explain the Monte Carlo method used for obtaining values of the 
magnetic scaling exponent. Simulations were performed on square lattices of size L x L 
with periodic boundary conditions, with L ranging from 4 to Lmax = 128 for g = 8, 64 
and Linax = 64 for q = 512. 

We employed the Wolff cluster algorithm [^]. The first part of the simulations was 
to determine the autocorrelation times r, which were found to increase with the lattice 
size and also with q. For the largest simulated lattices, we determined r as follows: 
88 ± 4 cluster updates for g = 8 and L = 128, 3000 ± 215 for g = 64 and L = 128, and 



g 


s* 


c 


c/log2(g) 


8 


6.5 (20) 


1.530 (1) 


0.5100 (3) 


64 


15.5 (20) 


3.050 (3) 


0.5083 (5) 


512 


32 (2) 


4.545 (10) 


0.5050 (11) 


4096 


65 (8) 


6.038 (24) 


0.5032 (20) 


32768 


135 (20) 


7.54 (3) 


0.5027 (20) 


262144 


250 (50) 


9.04 (3) 


0.502 (2) 



Table 2: Critical disorder strength s and central charge c, as functions of q. 

31000 ± 3000 for q = 512 and L = 64. This rapid increase of r with q explains why we 
simulate only up to L = 64 for the largest q. 

Next, we measure the magnetisation, defined for each disorder sample x by 

q{p) - 1 



m^ 



q-1 



(4.1) 



where p = max(A''i, N2, . . . , Nq)/L'^ and A^^- is the number of Potts spins taking the value 
a. Here (...) denotes the thermal average. Then the magnetisation m{L) is obtained by 
averaging over 10^ disorder configurations for q = 8, and 10^ configurations for g = 64 
and 512. For each disorder sample, 100 x r updates were dedicated to the thermalisation, 
and a further 100 x r to the magnetisation measurement. Error bars were computed from 



the disorder fluctuations (it can easily be checked [|10| that the contribution from thermal 
fluctuations is negligible), and the strength of the disorder was chosen as indicated in 
Table |. 

From a fit to m{L) ^ L~^% we obtain for the magnetic scaling dimension: 



Xi 



0.1535(10) for q = 8 
0.172(2) for g = 64 
0.180(3) for g = 512 



(4.2) 



We see that the magnetic exponent seems to saturate as we increase q. In view of 
the result (|1.1| ) for the central charge we expect the asymptotic behaviour should involve 
log(g) rather than q itself, and indeed the data are well fitted by 



xiiq) = a + b/ \og{q) 



(4.3) 



with a = 0.192(2) and b = -0.080(4). Thus, based on the form (^ we are led to 
propose the hmiting value ( |1.2|) of Xi given in the Introduction. 

5 Discussion 

It is useful to juxtapose our findings on the large-g behaviour of the critical line with the 
phase diagram proposed in Ref . |^ . In that work the disorder strength was parametrised 
through s = q"^ with w > 0, and the limit w —>■ oo was identified with classical percolation 
on top of the strong bonds. Actually it is easily seen from Eq. (p.2|) that directly at q = oo 
this percolation scenario holds true whenever w > 1/2, and assuming that the line of 
critical fixed points is described by a monotonic function w^:{q) it can thus be confined 
to the region w < 1/2. With this slight reinterpretation, Ref. §] argues that at q = oo 
the critical point is located in the limit w -^ 1/2. Indeed, since for q = oo any initial 
w <^ 1/2 will be driven to larger values due to mapping to the random field Ising model, 
this is nothing but the usual assumption of "no intervening fixed points" . 

However, this seems at odds with the results of Table ^ where we found that for g S> 4 
the critical line, when measured in terms of w, saturates at w = 0.31 ± 0.02. Unless our 
numerical method is flawed by some gross systematic error, it is thus a priori difficult to 
see how this can be reconciled with the above result of w*(g = oo) = 1/2. A possible 
explanation is that the limits q ^ oo and w -^ 1/2 are highly non-commuting. This is 
witnessed by the jump in the central charge, which in the percolation limit {w = oo and 

g < oo) reads [^ 

5\/3 
Cperc = -^Hq) - 0.47769 log2(g), (5.1) 

to be contrasted with our numerical result ( |1 . 1| ) . 
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